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generalize this argument to covariant derivatives with torsion. We find that some 
components of the torsion field can be identified with the dilaton and the i?-field 
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in which the equation of motion of string theory is exactly satisfied. 
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1 Introduction 



String theory, as a candidate of unified theory of particle physics, is expected to provide 
knowledge of physics governing high-energy regime and possibly early time of the universe 
where the quantum effects of gravity become important. It is of crucial importance to 
have descriptions of string theory which give insight beyond perturbative level. Matrix 
models give efficient ways of approaching nonperturbative nature of string theory [H [2] . 
In matrix models, how to describe the spacetime geometry has been one of key questions. 

In the original proposal of matrix models, matrices have been considered as spacetime 
coordinates. In this description, the gravitational interaction appears as a loop effect, 
but the origin of the general covariance is unclear. In the recent paper [3], another 
interpretation of matrix models was proposed by some of the authors. It has been shown 
that a covariant derivative on any d dimensional spacetime M can be expressed by a set 
of d operators that acts on the space of functions on the principal Spin{d) bundle over 
M. The Einstein equation follows from the equation of motion of the matrix model if 
we identify matrices with the covariant derivative in this way. Symmetries under local 
Lorentz transformation and diffeomorphism are manifestly realized as a part of the U{N) 
symmetry. On the other hand, a general matrix should be regarded as a power series 
in the covariant derivative and Lorentz generator that contains infinitely many fields as 
coefficients. To understand their dynamics, it is important to investigate the degrees of 
freedom which appear around a classical solution. In this paper, we study the fields that 
appear as the coefficients of terms linear in the covariant derivative and Lorentz generator, 
that is, the vielbein and the torsioijf]. We discuss a possible interpretation of these fields 
in the framework of string theory. 

The paper is organized as follows. In the next section, we review the new interpretation 
of matrix models [5]. In particular, we explain how the covariant derivative on a curved 
space can be expressed in terms of matrices. In section [31 we analyze the dynamics 
of the vielbein and the torsion. They form a closed subset of equations of motion and 
are expected to have consistent solutions without introducing higher spin fields. By 
identifying the trace and the antisymmetric parts of the torsion to be the field strengths 
of the dilaton and the antisymmetric 5-field, respectively, we find solutions which resemble 
those in string theory. However, it turns out that ghost like degrees of freedom can also 
propagate in general. We might need to impose a constraint on the expansion of matrices. 
In section HI we study modified models which have mass and cubic terms. In these cases, 
there are classical solutions with constant torsion. If we interpret the antisymmetric part 
of the torsion to be the field strength of the 5-field, then these solutions satisfy the 
Einstein equation with a constant background of field strength of B. This supports the 
idea of embedding the -B-field into the torsion that is introduced in section [31 Section 
[5lis devoted to conclusions and discussions on future directions. Some of the details are 
contained in the appendices. 

^ Higher spin fields appear from coefficients of fiigher powers of covariant derivatives, and tlieir gauge 
symmetries are included in the U{N) symmetry [3] . Free equations of motion are also correctly reproduced 
from that of the matrix model [J] . 
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2 New Interpretation of Matrix Model 



In this section, we briefly summarize the results obtained in [3]. We consider the 
large- reduced model of d-dimensional U (N) Yang-Mills theory, 



where A^ are N x N Hermitian matrices, and a and b run from 1 to d. (Supersymmetric 
version of this model with d = 10 is known as IIB matrix model [2].) This model has the 

U{N) gauge symmetry 

SAa = t[A,Aa], (2) 

where A is a Hermitian matrix. 

In [3] , a new interpretation is proposed in which the matrices Aa represent differential 
operators on a Riemannian manifold. In this interpretation, the symmetries of general 
relativity are contained naturally in the U (N) symmetry, and the Einstein equation follows 
from the equation of motion of the matrix model. 

To begin with, let us see how the covariant derivatives are described in terms of 
matrices. An important aspect of the covariant derivative is that it maps a rank-n tensor 
to a rank-(?7, -|- 1) tensor. For example, is not the a-th component]^ of V acting on 

the 6-th component of V, but is the (a, 6)-component of a rank-two tensor W; indeed, 
in the formula 



V,H = e/(9^H + ^^6^K), (3) 

all of the components Vc contribute to the covariant derivative. Here, e^^ and u^^'^ are 
the vielbein and the spin connection, respectively. 

Therefore, in order to express the covariant derivative in terms of matrices, the space on 
which matrices act must contain the tensors of all ranks. It was shown that the covariant 
derivative can be interpreted as a set of matrices (or equivalently, endomorphisms) acting 
on functions on the principal Spin{d) bundle over the Riemannian manifold M^. The 
space of functions on the group manifold Spin{d) forms the regular representation which 
contains tensors of all ranks. Then, the covariant derivative Va is mapped to a set of 
matrices V(a) by 

V(,) = R^a)\g-')^b 

= R^,)\g-')e,^{x) {d, + a;/"(x)ad) , (4) 

where g is the coordinate of the fiber Spin{d), R{^a)'{g~'^) is the vector representation of 
Spin{d) and Oat is the Lorentz generator of the Spin{d) group. The operator V(a) acts 

^ a,b, - ■ ■ represent the Lorentz indices. We use ^,v, - ■ ■ to express the Einstein indices. 
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as endomorphism on the space of functions on the principal Spin{d) bundle. Products of 
V(a) are related to those of Vq as follows: 

On the left hand side, each of the indices (a), (6), (c) is nothing more than a label and there 
is no mixing of them in the product. Thus the product can be represented as the product 
of matrices V(a), V(;,), V(c), ■ ■ ■. On the right hand side, product of Vq/, V;,/, Vc', ■ ■ • mixes 
the indices b',c',- ■ ■. 

We regard that the matrices Aa act on the space of functions on the principal Smn{d) 
bundle over M^. Then, Aa can be expanded as series in powers of V(a) and Oat, III 

Aa = a^a){x,g) + ^{6(^a)^'''> +b(^a)^'\x,g),V^b)} + ^{Cia)'"'{x,g),Obc} 
•2 '2 

+ '-{dia)^^'^^'^\x, ^7), V(,) V(rf)} + ^{{e(,)^'=W(x, g),Obc, }, V(,)} + ■ ■ ■ . (6) 

If we interpret the coefficients appearing in this expansion to be local fields, then the 
diffeomorphism, the local Lorentz transformation and higher-spin gauge transformations 
are contained naturally in the U{N) symmetry of the matrix model [31 Hj. For example, 
the diffeomorphism is generated by A = iX^^^ (x) V(a). (We will show the explicit expression 
for the diffeomorphism and the local Lorentz transformation in section [3l ) 
The equation of motion for the action ([1]) becomes 

[A'^,[Aa,Ab]] = 0. (7) 

As an ansatz, it is natural to require Aa to be the differential operators of first-order, 
because their commutator is again of first-order and hence they form a closed algebra. 
Here, we consider the simplest possibility among them, 

Aa = lV^a)- (8) 

(We consider a general ansatz in the next section.) Then, ([7]) becomes 

[V^ [V„ V,]] = {V^Rab'^)0,, - Rb'V, = 0, (9) 

where we have converted the indices (a) and (6) to the indices a and b using the equation 
dSD. Because WRah"'^ = follows from Rab = and the Bianchi identity V^^Rab""^^ = 0, the 
equation ([9]) reduces to Rab = and thus provides the Einstein equation. Furthermore, 
if we regard the fields appearing in ([6]) to be fiuctuations about the classical solution 
of the form ([8]), they can be interpreted as local fields propagating on the curved space 
Md which is represented by the covariant derivative V(a); for example, aa{x) satisfies 
the Maxwell equation [5] and da'^'^ satisfies the free equations of motion for higher-spin 
fields [1]. Therefore, we can interpret that the classical solution derived above describes 
the Ricci-fiat background spacetime. This is a new mechanism of the emergence of the 
spacetime in the matrix model. 

^ The coefficients a^a)ix, g),b^^J<^{x, g), ■ ■ ■ are real-valued. We introduce the anticommutator { , } in 
order to make Aa Hermitian. Hereafter, we omit the anticommutator for simplicity. 



4 



3 Torsion as Massless Fields 



In this section, we consider the general form of the first-order differential operator, 
which is given by 

A, = iD^,-) = tR^a)\g^')D, = iR(a)\g~') (e6^(x)V^ + S,''\x)0,d) , (10) 

where and Sa'"^{x) are the torsionless covariant derivative and the contorsion, re- 
spectivelylf] The differential operator Da can be regarded as a covariant derivative with 
torsionjj As discussed in the previous section, because differential operators of first-order 
form a closed algebra, it is expected that the equation of motion has consistent solutions 
by setting all the fields in ([6]) to be zero except for the fields appearing in the coefficients 
of ffTOj) . Note that it is possible to consider (yf-dependent fields eh^{x,g) and Sb"^{x,g) in 
general. In this paper, we restrict our attention to the simplest case where and Sh''^ 
has no (^-dependence. For details concerning this point, see [3]. 

Before analyzing the equation of motion, it is helpful to study the gauge symmetries. 
The action ([1]) is invariant under the U{N) transformation 5Aa = i[A, Aa], where A is an 
N X N Hermitian matrix. If we restrict A to be the first-order differential operator. 



A = ^A^(x)V^ + zA"''(x)a6, (11) 
then the local fields appearing in (ITOl) transform as 

5e/(x) = -Xa\x)ef,^ix) + V,A^(x), (12) 

5V(a^) = V^A^^(x) + X''{x)R,,''{x), (13) 

SS^ = -X^{x)V^Sa''%x) - X/S,'%x) - X\S/'{x) - X'dSa'^x). (14) 



Therefore, zA'^(x)V^ and iX°'''{x)Oab generate the diffeomorphism and the local Lorentz 
transformation, respectively, and we can use them to remove unphysical modes of Ca'^ix). 
Then, no component of the contorsion Sj"^{x) can be gauged away. 

Now, let us derive the equation of motion. The commutator of Da gives 

[Da, Db] = Uab'^O^d + Tab'D,, (15) 

where Tabc = Sate — Sbac is the torsion and TZat'^'^ is the Riemann tensor with the contri- 
bution from the torsion, which is related to the ordinary Riemann tensor Rab'^'^ by 

cd r) cd I v7 C v7 ccd q c q ed , q c o ed f ^ 

K-ab — J^ab + y aJb — V b^a — Oae Ob + Obe Oa ■ [i-O) 

^ Although wc can introduce cocfBcicnts of the covariant derivative a.s Aa = i (a(a)''''^ V(b) + S(^a)'"^Obc) , 
such coefficients can be absorbed by the foUowing field redefinition: e'^^-j^ = a(a)^''-'e(b)'', S'^^-j^^Obc = 

SabcObc + {^{a)^'''^^ [b) ~ ^(a)) ' H^rc Vjj is the torsionless covariant derivative which is constructed from 

■ 

^ The torsion gravity has been studied extensively as a generalization of general relativity, and it is 
interesting to clarify the relation of it to our approach. For reviews of the torsion gravity, see [SJ [5] . 
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The equation of motion is given by 



= 0, 



(17) 



where IZ, 



. Therefore, we obtain the following two equations: 



—TZbc + D"'Tabc + Tab'^T"'dc — 0, 



(18) 
(19) 



Note that f|T8l) has not only the symmetric part, which can be regarded as the Einstein 
equation, but also the antisymmetric part. These equations are solved in the following 
subsections. We first analyze them at the linearized level in the next subsection. We find 
solutions in which torsion fields can be interpreted as the dilaton and the antisymmetric 
5-field in string theory. In subsection 13.21 we show that the equations have solutions to 
all orders in a weak field expansion. 

3.1 Linear Approximation 

In this subsection, we solve the equations of motion f|T8l) and f|T9|) at the linearized 
level. The contorsion Sj^'^ can be decomposed to irreducible tensors as 



where Hahc is antisymmetric and Labc satisfies the relations L\^ahc] = and L"ac = 0. Using 
these irreducible tensors, we can linearize (fTS]l and ffT^ about the fiat space, and we 
obtain 



Note that here da represents 5a^d^. The expressions (1211) and fl22l) are the symmetric 
and antisymmetric part of ffTSjl . respectively. As we have discussed below ( IT^ . K, Habc 
and Labc are not gauge fields. Therefore, if we regard them as fundamental fields, there 
appear negative-norm modes0 Then, it is natural to regard them as field strengths. As 
we will see below, this assumption leads to equations of motion of massless fields in string 

^"^By negative modes we mean fields with odd number of time like index 0, although we have not started 
from any field theory action. It is hard to expect such degrees of freedom to posses positive norm when 
the equation of motion eould be derived from a field theory action. 



Sabc — Vh5ac — Vc5ab + Labc + Habc, 



(20) 




Rab = (3 - d)d(^aVb) - 2(9 ■ V)5ab + 29^L(,,), 
(l-rf)9[,H]+3ra,6e = 0, 



(21) 
(22) 
(23) 
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theory at least in the hnearized level. For this reason, we consider the following form of 
solutions: 



Va = da(p, (24) 
Habc = d[aBbc], (25) 
Labc = 0. (26) 

Then, the equations of motion reduce to 

Rab = (3 - d)dad,(l), (27) 
(9^0 = 0, (28) 
d^Habc = 0. (29) 

Here, the equation fl2I?l) is obtained directly from fl22]) . The equation (1251) is obtained by 
substituting fl^Tl) to the Bianchi identity V^-R = 2\/°'Rab- Using the equation fl25]) . we 
obtain the equation (127|) . Using the Bianchi identities V^°'Rab^'^ = and d[aHbcd\ = 0, we 
can see that (!23|) gives no other equation. 

These equations are the same as the linearized equations of motion for the system of 
the dilaton and the antisymmetric S-field coupled to gravity, which arise naturally from 
the string theory^ Note that the conditions ([MD, (125]) and are modified if we include 
the nonlinear terms. We discuss this point in the next subsection. 

Although we have considered solutions which satisfy the conditions fl2^ . fl25l) and 
f l26|) . as we will see in the next subsection, and in appendix |X1 there are many solutions 
which do not satisfy these conditions. Such solutions are undesirable, because negative 
norm states can not be removed by the gauge symmetry. We discuss this point further in 
section [51 

3.2 Weak Field Expansion 

In this subsection, we introduce a systematic weak field expansion. We find that for 
any solution of the linearized equations we can determine higher order corrections order 
by order. In particular, the solution we discussed in the previous subsection, which looks 
like the system of the dilaton and the S-field coupled to gravity, can be lifted to all orders. 

In order to perform the weak field expansion, it is convenient to redefine the fields as 

Da = ea^d, + S'a'^Obc, (30) 

where the quantity S'a^'^ is a sum of the contorsion and the spin connection, 

ST = Sa'" + ea'-UJ^"'. (31) 

By expanding e^^ about the fiat space as 

Ca^ = 6a^ + /l/, (32) 
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An earlier attempt to embed the i?- field in torsion can be found in [7]. 
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the equations of motion can be written as 

[D\ [Da, D,]] = {d'^U^ - j/) + {dPF^,''' - J,'") = 0, (33) 

where 

fab^ = dah^^ - d,Kr Fab"" = daS'C" - 8,3'^'^ (34) 

and ia^ and J^'^'^ are the polynomials of ha^ and 5'^'^'^. Note that all the indices a, 6, /i, z/, ■ ■ ■ 
are those in the flat space. In particular, the derivative da does not stand for Ca^d^ but 
Sa'^d^. The quantity Ja'^'^ has only the terms of degree two and three in ha^ and 5'^'^'^, 
while the quantity jb'^ has also the linear term; 

jb^ = -^d'^Sab^ + d^S'ba^ - ^dbS'a^^ + tcrms of degree 2 and 3. (35) 

Because 9^ and Ocd are independent operators, the equation (13^ is equivalent to the 
following equations: 

( [D^ [Da, Db]]\sr = d'^fab' - 3b' = 0, (36) 
( [D^ [Da, A]] la)'" = d^'Fab"' - Jb"' = 0. (37) 

These expressions look similar to Yang-Mills theory except that there are extra indices 
fi, c and d. 

Let us solve equations (15^ and (P7|) order by order in a weak field expansion. We 
expand /i^^ and S'a^^ as 

/i/ = + h^2)^. ^ . . . ^ ^/6c ^ gKl)bc ^ g,{2)bc ^ . . . ^ ^gg) 

where h^f^'^ and S'J^^^'"^ represent quantities of i-th order. The equations of motion at the 
z-th order are 

[Da, Db]f\Y = - {h('\ ■ ■ ■ , h^'-'\S'^'\ ■ ■ ■ S'^'-'\ S'^'^) = 0, (39) 

( [Z^^ [Da, Db]f\j' = d^F^"' - Jt' (/^W, ■ ■ ■ , h^^-'\S'^'\ ■ ■ ■ , ^'(-1)) = 0. (40) 



Note that jj^^^ contains the contribution from S'J^^'^^'^, while Ja^^^^'^ does not. 
At the lowest order (140|) becomes 

d'^Fii'"' = 0, (41) 

because Jj^^^^'^ = 0. This is nothing but the Maxwell equation with the extra indices c and 
d, and can be solved for S'J^^^'^'^. By substituting this solution S'j^^^'^'^ to fl39l) . we obtain at 
the lowest order 

d^fil^'=A'^'{S'^'^). (42) 



We can show that the divergence of jl^^^ vanishes automatically. Indeed, using fl35|) and 
[Tl). we have 



b-Wl^ f q'{l) 



ab 



0. 



(43) 



Therefore, (142)) has the same structure as the Maxwell equation with a conserved source 
except for an extra index /i, and has a solution for h^^^^^. 

In this way, we have solved the equation of motion at the lowest order. We next show 
that the higher order terms can be determined order by order. Suppose we have solved 
the equation of motion to the (n — l)-th order: 



b\r'=0 (z = l,---,n-l). (44) 



Under this assumption, we show that the equation of motion at the n-th order, 

[D^[D„A]]^"^=0, (45) 

or equivalently 



in) 



a 

in) 



ab 



O 



cd 



Jb 



ga-pin)cd _ jin)cd ^^(1)^ ^ h^"-^\ S'^'\ ■ ■ ■ , S'^^-^)) = 0, (47) 

also has a solution for h^^^^ and S'S"^^^^. What we need is to show 

0, 



. . . ^ /,(«-!) ^ ^'(1)^ . . . s'(^-^\ ^'W) = 0, (46) 



b 

b ■ 



d'j'b 



0, 



(48) 
(49) 



because if these relations hold, then fl46l) and fl47l) can be solved as in the case of fl42|) . 
First let us show (148|) . From (!34|) . d'^d^Pj^^^'^'^ vanish identically. Therefore, we have 



b j{n)cd 



{n)cd 



ab 



-J, 



(n)cd 



By definition, the right hand side of this equation is equal to 



d'([D-,[D,,Di, 



l(n) 



O 



cd 



(50) 



(51) 



Noticing that 



in) 



+d'([D'',[Da,Db]] 



in) 



O 



cd 



^in) 



on 



(52) 
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we obtain 



d'{[D^, [D^,D,]]^\y' = (- [d\ [D^, [D,,D,]] 



cd 



Then, using the equations of motion in the lower orders 



o 

we have 



d\[D^,[Da,D,]] 



(n) 











of- 


i- 



in) 



O 



cd 



(53) 



(54) 



Finally, from the identity D^, [Da-, -Df,]] = 0, which is easily derived from the Jacobi 
identity, the right hand side of (15^ is zero. This completes the proof of (l48l) . and hence 
we can find a solution of (147!) for S'J^'^^^'^. 

Next we show (HOl) . In this case we must use the equation of motion for S'^'*-"-""^, in 
addition to ones for /i^ and S'^^\i = l,---,?7, — 1), because j^"^^ depends on 5"*-"-* as can 
be seen from fl5^ . As before, using flM|) and flS2]) we have 



Jb 



i(«) 



cd 



OH 



(55) 



Then, by using (l47l) and the identity D'', [D°', [Da, -D^]] = the last expression becomes 
zero, and hence we have (l49!l . 

In this way, we can construct the solution order by order. Note that the resemblance to 
Yang- Mills theory is crucial for the existence of the solution. In appendix [XJ we discuss 
how many physical degrees of freedom remain after taking the gauge symmetries into 
account. 



3.3 Further Analysis for 

In section 13. ![ we showed that at the linearized level the equations of motion ffTSjl and 
([T9|) have a solution composed of the scalar and the rank-two tensor Bab- In section 
13.21 we showed that any solution at the linearized level can be lifted to a solution of full 
orders in the weak field expansion. Although the degrees of freedom other than (p and Bab 
are not desirable as propagating fields, they are necessary to make full-order solutions. In 
order to clarify this point, we discuss the difficulty of the following ansatz: 

K = V,0, Habc = 0, Labc = 0. (56) 

We show that (IS^ gives a rather strong condition on so that cannot be regarded as a 
scalar field. Note that we impose this condition to all orders. In this ansatz, the Einstein 
equation f|T8l) becomes 

Rab = (3 - d) Va V60 + (2 - rf) ( Va0) ( Vfe^) - 2( V ■ 4>)6ab + {2d - 3) ( V^)'^^^. (57) 
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Substituting it to the Bianchi identity 

2V"i?ab - VbR = 0, (58) 

we obtain 

da{id- 1)V20 + {-d^ + d + 3)(V0)2} 

+ (9,0) {{2d - 8)(V'0) + 2(3 -d){d- 1)(V0)'} = 0. (59) 

This condition is too stringent for to be a scalar field. Indeed, in the weak field expan- 
sion, a generic solution of the first order cannot be lifted to all orders. 
At the first order (l59l) reduces to 

920(^) = 0, (60) 

and hence (p^^^ is a harmonic function. Then, we determine e^^^^ from (157|) . Using (l60l) . 

becomes a quantity of the second order, and hence the second term in the 1. h. s. of 
(l59ll is a quantity of the third order. Therefore, at the second order (l59l) becomes 

{{d - 1) + i-d'' + d + 3)(V0)2}^'^ = c, (61) 

where c is a constant. Because (V^^)*-^^ can be rewritten as 

(V20)(2) =9^0(2) + (V20«)(2), (62) 

we can solve (j)^'^^ in terms of 0'-"'^^ and e^^-*'^, and then we can determine e'^^^ from (1571) . 
At the third order, fl59l) becomes 



= da {{d - 1)9^0(3) + (polynomial of 0^^^ 0^^^ e^^^ e^^))} 

+ (9,0«) {{2d - 8)VV + 2(3 -d){d- l){d^Yf^ 
= da {{d - 1)9^0(3) + (polynomial of (p'-'K e^^^, e^^))} 

+ (9,0W) \{2d - 8) ^':_^ ~^ + 2(3 - d){d - 1)| (90(^))^ 



da<P^'\ (63) 



2d -8 

where we have used ( l6Ti) at the second equality. In order for 0^^-' to exist, the second term 
of the last expression must satisfy the integrability condition, that is, (9,0*^^''^ ((90*^^''^ 
must be total derivative. However, this is not the case for generic 0*^^^ Therefore, we 
cannot set Habc = and Labc = ; even if we impose Hj^l = L^^J^ = 0, Habc and Late 
appears in higher orders as nonlocal quantities obtained from 0. 
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4 Matrix Model with Additional Terms 



In this section, we study matrix models with a mass term and a cubic term, which 
were introduced in the context of D-brane action in background fields [H [9] or as toy 
models to obtain nontrivial vacua [101 [H]. We analyze the effects of these terms in our 
interpretation. In section 14. 1^ we consider matrix models with a mass term. In the 
paper |3], it was shown that Einstein manifolds are obtained as classical solutions of 
these models, and for such solutions the coefficient of the mass term plays the role of the 
cosmological constant. If we introduce a torsion, we find another type of solutions that 
have constant background of field strength of B. In this case, the coefficient of the mass 
term is given by the field strength. In section 14. 2[ we consider a three-dimensional matrix 
model with a cubic term. We find that this model also has the second type of classical 
solution in which the coefficient of the cubic term can be interpreted as the field strength. 
In both cases, the antisymmetric part of the torsion equals the field strength of B. This 
result supports the interpretation that the 5-field is embedded in the torsion. 

4.1 Matrix Model with Mass Term 

In this subsection, we consider the matrix model with a mass term 

S = -^Tr[Aa, A'] - ^S^bTr (AM^) , (64) 

where a is a parameter with the dimension of mass. We show that this model has two 
types of solutions without and with torsion. The former consists of the Einstein manifolds 
and the latter does group manifolds with constant field strength of B. 

First, let us consider the case where the torsion is absent [3]. For an ansatz Aa = ^V(a), 
the equation of motion 

[A^[A„,A]]-a'^ = (65) 

becomes 

[V^ [v., Vb]] + a'Vb = - [Rt' - a^6t') + {V''Rab'''')0,d = 0, (66) 

and this is satisfied if the spacetime is an Einstein manifold, 

R^^ = a^Sb". (67) 

The quantity plays the role of a cosmological constant. 

We next consider the case in which the torsion is present. For an ansatz Aa = iDf^a)-, 
where Da is a covariant derivative with a torsion, the equation of motion flB^ again 
becomes 

[D\ [Da, Db\] + a^Db = 0. (68) 
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If the operator Da is proportional to the standard left-invariant derivative on a Lie group 
G, they satisfy the commutation relation of the Lie algebra, 

[Da,D,]=afatcD,, (69) 

where fabc is the structure constant of G. Then the equation of motion is satisfied if 
fabc is normalized as 

facdfb"'^ = Sab- (70) 

As shown in appendix [Bl we can interpret Da as the covariant derivative on G with 
torsion, 

b ^ 
-Da = Va + Sa '^Obc, Sabc = Habc = -^fabc- (71) 

Here, the vielbein is defined by 

ta = e/ ■ tag~'d^g, (72) 

where g{x) is an element of G, x is a coordinate on G and ta is a generator of G, and we 
have 

r? cd " r rede tt rjcde 

J^ab — ^JabeJ — ^abeJ^ , ['O) 
Va^bed = VaHbcd = 0. (74) 

The torsion (ITT!) can be regarded as a constant field strength of Bab, 

Habc = ^[aBbc]- (75) 

Indeed, from (173!) and (1741) we obtain the equations 

Ra' = HaefH''^, (76) 

V^i/afe, = 0, (77) 

which are nothing but the equations of motion for the metric tensor and the S-field 0. 
Note that we could set Va = and Labc = in this case while it is not the case when the 
mass term is absent. 



4.2 Matrix Model with Cubic Term 

In this subsection, we show that the three-dimensional matrix model with a cubic term 
also has the group manifold SU (2) as a classical solution. 

This set of equations can also be regarded as equations for the graviton, the _B-field and a constant 
dilaton with a cosmological constant A = -^cP- 
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Let us consider the action 

S = -^Tr[Aa, A'] + ^eabcTr[A-, A']A'^, (7J 

where eabc is the structure constant of SU{2). For an ansatz Aa = iD{a)i the equation of 
motion is 

D(a), [^(a), ^(6)]] - aefccd[£'(c), ^(d)] = 0. (79) 

Because the structure constant eabc is an invariant tensor of the Lorentz group Spin{3) = 
SU{2), we can rewrite (1791) as 

[Da, [Da, A]] - (yebcd[Dc, D,] = 0. (80) 

Therefore, if operators Da satisfy the commutation relation 

[Da, Dh] = aeabcD,, (81) 

then Aa = iD(^a) is a classical solution. As in the case of the previous subsection, this 
solution describes the group manifold SU (2) ~ with a constant field strength of B. 

We now comment on dynamical properties of classical solutions. The solution Aa = 
iD(^a) can be expressed by a Lie algebra corresponding to the group manifold. Therefore 
we can discuss dynamical properties of this background using techniques of the noncom- 
mutative geometry. On the other hand, the solution without the flux does not form a Lie 
algebra, and hence the analysis is rather difficult. As a concrete example, we discuss a 
stability of in appendix [Cl 



5 Conclusions and Discussions 

In the previous paper, we have introduced a new interpretation of the matrix model 
where matrices act on the space of functions on the principal Spin{d) bundle over a Rie- 
mannian manifold. In this paper, we studied the equation of motion for the vielbein and 
torsion in the matrix model. In section [3l we interpreted the trace and antisymmetric part 
of the torsion to be field strengths of the dilaton and the -B-field, respectively. We found 
a solution in which only the graviton, dilaton and -B-field propagate. In general, however, 
all components of the torsion, including the negative-norm modes, can propagate, and 
there is no gauge symmetry to remove them. In spite of this difficulty, there are several 
possibilities that these unnecessary fields disappear from the spectrum. First, recall that 
higher spin fields couple to the torsion fields in a complicated manner, although we have 
simply set them to zero. They may give the large mass for the unnecessary degrees of 
freedom through the interactions. Also there is a possibility that the mass terms of un- 
necessary fields arise through quantum corrections. If the unnecessary fields disappear 
through such mechanisms, then our system can reproduce the massless bosonic sector 
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of superstring theory. In section HJ we studied the matrix model with mass and cubic 
terms, and found solutions with a constant torsion. In this case, the antisymmetric part 
of the torsion can be interpreted as the constant background of field strength of B. This 
result adds another ground to regard the antisymmetric part of the torsion to be the field 
strength of the 5-field. 

If unnecessary fields do not disappear from the spectrum dynamically, we may need to 
impose certain constraints to eliminate them. Although to find a good constraint seems 
difficult, it might become easier if we first extend the matrix space. One possible way 
to extend the matrix space is to replace usual manifolds with supermanifolds [13]. This 
corresponds to consider supermatrices instead of ordinary matrices. Due to this extension, 
we can implement the local supersymmetry as a subset of the unitary symmetry. In this 
case, a supermultiplet including the gravity is embedded successfully and the torsion 
constraint removes unnecessary components of the torsion. In particular, in the case of 
four dimensions it was shown that the equation of motion of the supermatrix model is 
consistent with that of A/" = 1 supergravity. This result suggests that a good constraint 
such as the torsion constraint can be found by extending the matrix space. Although 
in ten dimensions the superspace formulation is rather complicated and it is not clear 
whether the torsion constraint fits well to our scheme, the research along this line would 
be helpful to understand how the spectrum of string theory is derived from the matrix 
model. 
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A Physical Degrees of Freedom 

In this section, we consider the physical degrees of freedom of the torsion. We can 
impose the gauge fixing condition 

d^S'^'"' = 0, d^'ha" = (82) 

at all orders in the weak field expansion. This condition removes the components along 
the light-cone direction, h^^ and S'_^'^'^. Furthermore, because the equation of motion 
f l36|) and fl37j) can be regarded as the equation of motion for Yang-Mills theory with extra 
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indices /i, c and d, we can use the residual symmetries to remove h^^ and S'_^'^. Therefore, 
the physical degrees of freedom are 

Si'". (83) 

Here, i represents the transverse direction, i = 2,3, ■ ■ ■ ,d — 1. In terms of irreducible 
tensors with respect to the transverse SO{d — 2), they can be decomposed as follows: 

Q' Ox.. Q' lA.. 

QfOl Q/iO QUI Qlij Qf jkk) n/ 1 "^{kk) 

) ) ) ' ^(ij) ^ _ 2 ' '-'{ij) d — 2 ' 

^[ij] > '^[ij] > -^fo-fc]. - ^_ 2 ~ ^[ijk] y^^) 

Unlike in the Yang-Mills theory, fields with indices and 1 are not completely removed. 
Quantities including as indices, hi^,Sf^ and S^'^^, are negative- norm modes. 



B Differential Calculus on Lie Group 

In this section, we provide basic results concerning the differential calculus on a Lie 
group G, which are necessary in section HI For further details and other applications to 
physics, see [H] for example. 

Let ta be the generator of a Lie group G which satisfies the relations 

[ta,h] = ifabctc- (85) 

The left-invariant derivative Da on G is defined by 

e''DJ{g) = f{g{l-ze^Q)-f{9), (86) 

where f{g) is a function on G. From (!85|) . we can easily see that Da satisfies the commu- 
tation relation 

[Da,D,]= fabcD,. (87) 

By introducing a local coordinate x. Da can be expressed as 

Da = Va^d^, (88) 

where 9^ = d/dx^ and Va^{x) (a = l,---,dimG) are vector fields on G. By taking 
f{g) = g, Va^ can be expressed as 

ta = Va^ ■ig~\x)d^g{x), (89) 

where g{x) is the element of G specified by the coordinate x. 
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Next we identify to a vielbein e^^. From this vielbein, we obtain the standard 
covariant derivative = e^^ (^(9^ + u^^^O^^, where u is the spin connection. From fl88l) . 
Va is related to Da as 

Da = ya- UJa'O^c- (90) 

Then, if we regard the second term of the r.h.s. of (190|) as contorsion, Da can be regarded 
as the covariant derivative with torsion, 

Da = Va + Sa^Oiic, Sabc = —^abc- (91) 

Combining it with (187|) . we have 

nab"" = 0, (92) 
Sabc = -UJabc = ^fabc (93) 

We can show the relation V ahcd = by using the Jacobi identity: 

"^^aSbcd = afbcd 

= Ga^dfj^fbcd + ^ab" feed + ^ac" fbed + ^ad" fbce 

2 ^fo-b feed ~l~ fae fbed ~l~ fad fbee) 

= 0. (94) 

By using (ITSl) . (fT6|) . (!92l) and (l94l) . the Riemann tensor can be expressed in terms of the 
structure constant (or equivalently, the torsion) as 

Tj ed T-7 Q ed I T-7 Q ed I c e c ed n e q ed 

J^ab — —^a^b + ^ b^a + ^ae ^b ~ ^be ^a 



^ f f c n ed r ere 

^ yjae Jb Jbe J a 

-Jabef""'. (95) 



C Group Manifolds in Matrix Model 

In this section, we discuss the covariant derivatives on group manifolds. We begin 
with the discussion for general group G and then discuss the case of SU (2) in detail. 



C.l Covariant Derivatives on Group Manifolds 

In this subsection, we consider the covariant derivative on a group manifold G. As 
shown in preceding sections, the left-invariant derivative can be identified with the covari- 
ant derivative Da with an appropriate field strength of B. The Hilbert space V on which 
Da acts is the set of functions on the principal S'pm((i)-bundle over G, where d = dimG. 
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As we can see from (1891) . the principal bundle can be covered by a single patch (i.e. the 
index a in (1891) is defined globally), and hence this bundle is the direct product of G and 
Spin{d). Therefore, we have 

V = C^{G X Spin{d)) = C^{G) ® G°°{Spm{d)). (96) 

From (1861) . Da acts on the first component G°°{G) as the left-invariant derivative, or 
equivalently as the regular representation. If we denote the regular representation of the 
generator ta by Tj^'^^\ Da can be expressed as 

Da = -iT^^'^s) ^ I (97) 

Using the standard procedure explained in section [2l the index a can be rewritten to that 
with a parenthesis: 

D^a) = Ria)\9'')D, = -iTi'-'^ ® Ria)\9-'). (98) 

A remark is in order here. Although Da satisfies the commutation relation of the Lie 
algebra, D(a) does not for G ^ SU{2) = Spin{3) because the structure constant fate is 
not invariant under the action of Spin{d). Indeed, the commutation relation of D(a) is 

= R(^a)''\g-')Ri,)'\g-')[Da',Dy] 
= R(^af'{9-')Rib)''i9-')fa'b'c'D,, 

= Riaf'{g-')R(b)''{g-')Ric/ig'')fa't.'c'D^,), (99) 

where g G SU{2), and the coefficient R(a)"'' {g~^)R{b)^' ig~^)R{cf' {g~^)fa'b'c' depends on g 
unless G = SU{2). In the case of G = SU{2), this coefficient does not depend on g, 

R[af\g-^)R[b)'\g-^)R[cf\g-^)ea'b'c' = ^abc. (100) 

and D(a) also forms a Lie algebra. Therefore, we can find the explicit form of D(a) in 
terms of matrices. 



C.2 Explicit Form of the Covariant Derivative on SU{2) Group 
Manifold 

In this subsection, we show the explicit form of the covariant derivative on group 
manifold SU{2) in terms of matrices. Let D'^ be the covariant derivative with torsion 
which is obtained from Da, 

Da = Va + ^eabcObc, (101) 

by fiipping the sign of the torsion: 

D'a = Va- leaMc- (102) 

18 



The basic observation is that Z)(a) and Z^l^-, commute with each other and that both of 
them satisfy the commutation relation of SU{2): 



Therefore, using an appropriate basis we can write and as 



-eabcD'^c)- 



D 



(a) 



Ijireg) 



(103) 



(104) 
(105) 



where J^'^'^s) jg ^^j^e regular representation of SU{2). The explicit form of J^^^s) jg given by 



j{reg) 



\ 



J[i/2] ® h 



J[l] (g) 121+1 



V 



(106) 



where J[/] stands for the spin-/ representation. Using (11041) . we can express the covariant 
derivative without torsion V(a) and the Lorentz generator Oab as 



ab 



(107) 
(108) 



In the latter part of this subsection, we discuss whether or not the background 
Aa = iD(^a) given in the previous subsection is stable in three-dimensional matrix models. 
Although the physical interpretation is different, this background is mathematically the 
same as multiple fuzzy sphere configuration arranged in a concentric pattern, whose dy- 
namical properties have been studied extensively. In the model defined by the action ( |78l) . 
this configuration has flat directions J^^ —>■ J^^ + c" ■ I21+1, where c" is a constant which 
represents the translation of each fuzzy sphere, and furthermore, if Ca takes a non-zero 
value, then tachyonic modes arise and the multiple fuzzy sphere configuration decays to 
a single fuzzy sphere which is described by the irreducible representation [T5j. There- 
fore, the background considered here is unstable in this case. This instability may be 
removed by introducing supersymmetry or by adding a mass term. For example, the 
three-dimensional supersymmetric matrix model with the cubic term is considered in [TU] 
and it was shown that multiple fuzzy sphere configuration is BPS if they are arranged in a 
concentric pattern. Therefore, the configuration fll04p would be protected from the quan- 
tum correction and be stable, while other configurations which are shifted with nonzero 
are not. In the case that the action has the mass term, it is easy to see that there is 
no flat direction. In [H] , three-dimensional matrix model with the cubic and the positive 
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mass term is considered at one-loop level, and it was shown that the multiple fuzzy sphere 
configuration becomes the true vacuum!^ The multiplicity of each irreducible represen- 
tation depends on the mass and the coefficient of the cubic term, and it is interesting to 
investigate whether the configuration fll04p can be a true vacuum. 
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